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Abstract 

Studying loop corrections to inflationary perturbations, with particular emphasis on infrared factors, 
is important to understand the consistency of the inflationary theory, its predictivity and to establish 
the existence of the slow-roll eternal inflation phenomena and its recently found volume bound. In 
this paper we show that ("-correlators are time-independent at large distances at all-loop level in 
single clock inflation. We write the n-th order correlators of ( as the time-integral of Green's 
functions times the correlators of local sources that are function of the lower order fluctuations. 
The Green's functions are such that only non- vanishing correlators of the sources at late times can 
lead to non-vanishing correlators for ( at long distances. When the sources are connected by high 
wavenumber modes, the correlator is peaked at short distances, and these diagrams cannot lead 
to a time-dependence by simple diff. invariance arguments. When the sources are connected by 
long wavenumber modes one can use similar arguments once the constancy of ( at lower orders was 
established. Therefore the conservation of C at a given order follows from the conservation of ( at 
the lower orders. Since at tree-level ( is constant, this implies constancy at all-loops by induction. 



1 Introduction and Summary 

This is the sequel to a series of recent papers [1, 2, 3] that we wrote studying loop corrections 
to inflationary observables, with the purpose of understanding possible large infrared effects that 
might arise when considering the theory of inflation at quantum level. In the latest of these [3] , we 
proved at one-loop level that in single clock inflation, that is inflationary theories where there is only 
one relevant degree of freedom, the correlation of the curvature perturbation ( is time-independent 



on distances much longer than the horizon. Studying loop effect is in general important for the 
consistency of the theory at quantum level. However, the time-independence of ( is very important 
also for the predictivity of inflation and for the existence of slow roll eternal inflation and its recently 
found universal volume bound [4, 5, 6] (see the introduction of [3] where we recently discussed in 
detail about these points). 

The purpose of this paper is to extend the former calculation to all-loop orders. We show 
the constancy of ( iteratively in the number of perturbations. At linear level it is well known 
that ( becomes time independent after horizon crossing. Then, assuming ( is constant at a given 
perturbative order, we show that ( is constant at the next perturbative order which provides a proof 
by induction. We achieve this in the following way. The non-linear nature of the equations forces a 
coupling between different modes. The nature of the Green's function for the ( fluctuations forces 
any interaction that happen in the past not to induce any time-dependence. Time dependence 
can only come from interactions that happened during the last Hubble time. Further, as it has 
been already shown in several ways (see for example [7, 8, 9, 10]), all interactions involving modes 
longer than the horizon are vanishingly small. This means that in the latest possible interaction at 
least two short wavelength fluctuations must be present, plus an additional arbitrary number of long 
wavelength fluctuations. This reduces the problem to computing correlation functions of local sources 
made of short-wavelength and long-wavelength fluctuations. The long wavelength fluctuations are 
lower order in the fluctuations, and so they can be assumed to be constant in time. However a 
constant-in-time long wavelength fluctuations is locally unobservable and therefore it cannot affect 
the short wavelength fluctuations. This means that the correlation of sources is dominated either by 
correlation of short modes at far away distances, that are too small to induce a time-dependence, or 
by correlations of derivatives of lower-order long-wavelength fluctuations, that decay in time. This 
implies that ( goes to zero at late times, proving the constancy of ( at all loop order in single clock 
inflation. In other words, ( goes to zero in single clock inflation as an operator ^. 

2 Local gauge for a long mode 

It will turn out to be useful for the proof of the constancy of ( to introduce the following set of 
coordinates, which was introduced for the first time in [11, 3] at linear level in the fluctuations. Here 
we generalize them to the non-linear level. If we consider a set of ( fluctuations that involve, at 
the time of interest, long fluctuations characterized by the highest wavenumber fciow already very 
outside of the horizon, we aim here to provide a change of coordinates valid in a local patch of 
volume V much smaller than and much larger than H~^, that takes us from the metric written 
in standard ("-gauge to a form that is locally of the form of a homogeneous anisotropic universe, up 
to corrections that go as {kf^^Vy^^. We wish to find a set of coordinates where this fact is manifest. 
Let us start in standard ("-gauge, defined by setting the scalar field to be uniform (50 = or the 
Goldstone boson tt = 0), and the metric to take the form 

ds^ = -N^dt^ + hij {dx' + N'dt) {dx^ + N^dt) , (1) 

^As described in [3], the non- linear dependence of ( on lower order scalar or tensor fluctuations is para- 
met rically different. This implies that the analysis can be done neglecting tensor modes, as we will do, though 
we expect that very similar arguments apply also when tensor modes are included. 
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in the ADM parametrization. The spatial metric hij takes the following form 

hij = 6ij e^^ . (2) 

The variables N and iV^ solve some constraint equations. For standard slow roll inflation (but the 
argument can be straightforwardly generalized to all models described by the EFT of Inflation [12]), 
they are the following 

MliVi [iV-i (^E'j + d'jE^i)] = , (3) 



where R^^"^ and V are taken with respect to the spatial metric hij. What makes the long wavelength 
limit in ("-gauge non-trivial is the fact that at linear order in Fourier space the solution for N'^ in the 
long wavelength limit reads 

N' ^ §^a(fc) , (4) 

^low 

which does not have a nice fciow limit for a generic time-dependence of (l- The same holds 
at higher order. If we go back to real space, this shows that at linear level N'^ oc x'^: N'^ becomes 
larger and larger at longer distances. The reason why this is not a problem is that N'^ is not really 
a physical observable: as we will show, what matters is rather djN^ . 

In order to show this explicitly, we will solve the constraint equations (3) around the origin, 
at leading order in {kf^^VY^'^. We can Taylor expand C^L^x^t) = CL(0,t) + 5^Clo^^ Because of 
rotational invariance, we can choose the gradient of Ql to point in the z direction. It can be checked 
that, non-linearly in C,l^ but at leading order in {k^^^VY/^ ^ we have 

Ar = i + ^ + o((fc,„^F)V3) , (5) 

where all the are evaluated at the origin and A^^^o are three space-independent functions. For 
later convenience, let us define 

^(^-) = 2^^Ta^'- 

As mentioned, the problematic term is A^^: it does not appear to go to zero as /ciow ^ 0, even 
after it has been acted upon by a spatial derivative. Notice however that this term diverges linearly 
in suggesting that the physically sensible quantity is rather djNi. Indeed, the structure of the 
constraint equations show that one can restrict oneself to solutions that diverge at most linearly 
in the spatial coordinated. This means that in order to find a reference frame where the universe 
looks locally manifestly like an homogeneous anisotropic universe, it is enough to find a coordinate 
transformation that set to zero locally Ni and djNi^ with the remaining fluctuations going to zero 
as fciow 0, after have been acted upon by a spatial derivative. We will work without assuming 
any specific time dependence for so that the following change of coordinates holds also when one 
cannot use the solution to the linear equations for as it happens for example when computing 
loop corrections to inflationary correlators [3]. 
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We are going to perform the following change of coordinates: 

x' = M'jit) + &{t) , with det[M] / , 



(7) 



which keeps the inflaton fluctuations zero. Under this change of coordinates, and expanding around 
the origin, the metric takes the form 

ds'^ = -N^dt^ + e2C(^(^'^)'^)5^^- [m^ dx^ + (^M'lX^ + C' + N'\^ + diN\ [M^mx"^ + C^)) dt 

X M^ndx''+(^Wnx'' + C^ + N^\^ + dnN^\^{M''pX^ . (8) 

To bring the metric in our desired form, we therefore need to solve the following two equations: 

& + diN\& + N\ = Q , (9) 



The generic solution reads 



11^ 



C\t) = - [u{t)-'] 



fdt'U{t'] 

Jt2 



(10) 



M(ir,= {V{t)-'-\\M{, , 



where 



(11) 

with T staying for time-ordering, ti^2^Ci^ Mlj representing integration constants, and dNl^ is the 
matrix such that [dNl^Yj = djN^\^. 

After substitution, the metric takes the form 



ds'^ = -N^dt^ + e2«^(^'*)'*)5,.M^M^'^ dx^dx"" 



(12) 



Since M contains only dN terms, the metric goes to zero as fciow — ^ 0, if acted upon by spatial 
derivatives, as we wished to show. The local spatial metric becomes 



space ^ ^2 2a(t) 1+En=l,...,5 ^I^,n W Sn,ij 



(13) 



where (l is the local ( fluctuation, that multiplies the identity matrix 1, and AL,n's are the local 
anisotropy coefficients, represented as coefficients of five independent symmetric traceless matrixes 

o 2 

Our choice to align the gradient of (^l with the z direction simplifies formulae remarkably. We 
have: 

/ 1 \ 



[dN\o]ij = SisSjsFiCL) . U = 



1 

Vo e^A^^^^(^-(^^)) J 



^ These are given by 



1 . . 2 _ 1 



(15) 



(14) 



Si 



5,U 



1, , 1, , 2^ ^ 



where (i, j) means symmetrization with respect to the i, j indexes. 



In the local frame the fluctuations, represented by a tilde, are given in terms of ^ by the following 
relation 



~ ^ 1 /•* 

CL{x,t) = CL{S{x),t) + - dt' 



H 2iJ + a(0,O 



a(o,t') , 



(16) 



/' 



g 2fl+a(o.«') ^ - 



All the remaining Al,^ vanish. For this reason, in the rest of the paper, we will call Al,5 as simply ^. 
Notice that under this change of coordinates only the long wavelength C,l is transformed. Short 
wavelength fluctuations ("5 transform as a scalar. The long wavelength metric takes locally the 
following form 

PL,M. = -N^dt" + g'l^'' dx'dx^ , (17) 
with the spatial metric defined in (17) and given here in matrix form 



/ g2(a(t)+|Ai(t)) 



fr%t) = a(tf 










g2(a(t)+|Ai(t)) 





\ 

g2(a(t)-|Ai(t)) J 



(18) 



and 



N{t) 



Cl (Ut)) 



(19) 



Notice that det[^L] = N^a^e^^^. The long wavelength mode has been locally reabsorbed in an ho- 
mogenous anisotropic background up to corrections that are manifestly subleading in (kf^^V)^/^ 
limit. Short wavelength fluctuations live in this very smooth universe. 



3 Proof of the conservation of ( to all loops 

Let us consider a long wavelength fluctuations (lix^t), defined as 

a(f , t)^ J d^'k e(A;iow - k) e''^-\j:it) , (20) 

where k\o^ is some comoving momentum. We are interested in proving the constancy of correlation 
functions of C,L{x^t) once all the comoving modes that are part of its spectrum are outside of the 
horizon. As we solve iteratively in the interactions the equation of motions for (^7,, we need to 
insert interactions at some specific times and spatial locations, and sum over these. In [1, 3], we 
have shown how solving iteratively in the interactions the equations of motion and then taking 

^ Later in the proof of the constancy of we will need to consider configurations where there is a generic 
lower-order long- wavelength configuration sourcing a long wavelength fluctuations. A priori, these might have 
different gradient directions. However, since we work at linear order in the gradients, rotational invariance 
implies they must have the same direction. Alternatively, one could imagine to perform a rotation to bring the 
grandient in the z-direction, obtain the local frame in this simple coordinates, and then perform the inverse 
rotation. In this way, one would obtain all the Xl^w In any event it is straightforward to use the generic Xl^u 
from the get-go. 
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correlation functions of the non-linear solutions is equivalent to computing loop corrections. Given 
a certain solution with the insertion of a given number of interactions, there will be one interaction 
that has been inserted last. Let us consider it. It is sufficient to work with the tilde variables of 
the local frame defined in Sec. 2. We can write the equation of motion for (^l a.nd the anisotropy 
coefficient in the following way 

{dt + m c!£\x, t) = {dt + m) (^^^ Cir^) = t) , (21) 

where 4^^'""^ and ^^^'""^ are some source that depend in a local way on (l^ A^, SNl^ and their 
derivatives, where Cl, 5Nl^ Al are respectively the time-shift and the anisotropy coefficient of the 
local space-time described in Sec. 2. Sl additionally depends in a more unconstrained way on the 
additional short wavelength fluctuations. Let us explain how we obtain the above equation: 

1. First of all, notice that the differential operator acting on (l on the left is very similar to 
the differential operator that acts on linear fluctuations in the standard ("-gauge, apart for a 
V^/a^ term which is missing. The gradient is indeed negligible because we are interested in 
proving the constancy of when fciow is already very outside of the horizon. Since, as we 
will see next, both the ("-Green's function and ("-Green's function of the long mode makes any 
solution that is sourced by a source in the past approach a constant, it is only the interactions 
that are happening at the latest possible time that could lead to a time-dependence. This 
means that in (21) we can restrict ourselves to times when fciow is outside of the horizon, and 
so the term V^Cl is negligible. 

2. The right hand side of (21) contains the source term Sl- As we said, this term is a local 
function of ("l, A^, SNl and their derivatives, and of the short wavelength fluctuations in 
a more unconstrained way. In general, each of the fluctuations in Sl will be either a free 
field, or the results of some interactions. Each interaction is associated to a vertex with 
some fluctuations, and so it is useful to express the perturbative series in terms of number of 
fluctuations involved, rather than in the number of vertices included. This is equivalent to the 
counting of loops ^. Let us call (^^"^ as the solution to the perturbative equations of order n 
in the fluctuations. If we wish to source the sum of the orders of all the (l, A^, SNl and 
the short fluctuations in Sl must be equal to n. In particular, notice that the fluctuations 
Cl: and 6Nl appearing in Sl are of lower order in the fluctuations. The structure of the 
equation suites itself therefore to an iterative solution. We will prove that the assumption that 
the lower order (l is constant implies that the next higher order (l is constant as well. Since 
at tree level (l is constant, the argument provides a proof to all orders ^. Notice immediately 
that assuming the lower order (l to be constant simplifies the structure of the equation, as if 
(l is constant, 5Nl must vanish and Xl must go to a constant value. 

^The same phenomenon happens in standard quantum field theory. For example, in the correction to the 
propagator at first order in diagrams with either one quartic or two cubic vertices contribute. In standard 
quantum field theory, indeed, fluctuations are of order h. 

^When ( is not constant at tree- level or alternatively when one cannot go to standard (-gauge as defined 
in (2), as for example in [13], our proof does not apply. 
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3. How can we guarantee that Sl is spatially local? Once a mode is much longer than the Horizon, 
the Green's function associated to the differential operator in (21) is local in space, as we will 
show next in (23). This means that a long wavelength mode at one location is sourced only 
by a source at the same location. This shows that the interactions of the long mode with the 
other modes can be effectively described in the local anisotropic frame introduced in Sec. 2, 
and we can use directly the variables with tildes. Notice that the constancy of Cl, ^L,n implies 
the constancy of (l and viceversa. In this frame it is manifest that the dependence of S'l on 
the long- wavelength fields is manifestly local. In more detail it follows from the fact that in the 
constraint equations for the short-wavelength SNs and iV^ there cannot be spatial derivatives 
acting on (l and as they can be taken as homogeneous. Corrections to this are explicitly 
suppressed by powers of k\o^/{aH). Notice that this is not the case for standard ("-gauge, 
where operators such as di/d'^ acting on long- wavelength fields, appear in the source. 

Since we are interested in the solution to (21) can be written directly: 

Ci"^(x,t)= / rfVP,,GL(x,xO;5f^'"\xO , (22) 
Jv 

where Gl{x^ x') is the ^L-Green's function for the long mode 

Gl{x,x') = 5('\x - x') GlM , GlM = ^ 3ff ' (23) 

in the slow roll approximation, and V^/ is the x'-differential operator that acts on ( in the vertex 
that gives rise to Sl- Since we can neglect gradients of the long mode, the Green's function is local 
in real space. We have 

CP (t) = I dt' A' GLit, t') 5f it') . (24) 

As promised, for t ^ t\ Gl tends to a constant, which means that in order for a time-dependence 
to be induced, we need to concentrate on t' ^ t. All of the above applies identically to the sourcing 
of Al after replacing (l with Xl apart from an irrelevant multiplicative constant in the definition of 
the Green's function, that is time-idependent at leading order in slow-roll. Slow-roll corrections are 
irrelevant for the purpose of proving the constancy of (l- 

Since generic n-point functions are formed by decomposing them as product of connected cor- 
relation functions involving a smaller number of fluctuations, we can concentrate on connected 
correlation functions, that are therefore simply related to connected correlation functions of Sl- 
Since the case of correlation functions involving also Xl follows virtually unaltered from the case of 
the correlation function of Cl? we quote explicit results only for (l. We concentrate first on the two 
point function. We summarize the following points: 

1. Some of the operators contributing to at least one of the must involve short modes, as 
otherwise all interactions vanish. 

2. We are interested in contributions at late times because the form of the Green's function 
makes (l go to a constant if the final interaction is in the distant past. 
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3. Any vertex connected to the x' by a short mode must be within a space-time distance of 
order Hubble at most. This is because short modes need to be oscillating and they become 
uncorrelated on larger distances or longer times ^ . 

Let us consider the correlator of S^^^'^\ and first consider taking expectation value only over the 
short modes. This expectation value should be thought of as taken on a diff-invariant regularized 
theory, for example in dimensional regularization, with the counter terms provided naturally by the 
operators present in the effective theory of single clock inflation [12] ^. The general structure takes 
the form: 

(25) 

where . . . represents A/,, 5Nl and higher time derivatives of all these functions evaluated at the same 
location. The first term accounts for cases where the vertices around xi and X2 are connected by high 
momentum lines (which there must be at least two by momentum conservation). Note that if there 
are short momentum lines connecting the vertices they have to be close both in space and in time. 
We call these Short-Connected {SC) diagrams and they reduce in the two point function at 1-loop 
to the cut-in-the-middle diagrams discussed in [1, 3]. Since at large distances these modes become 
uncorrelated, the contribution is peaked at the same spatial location by a simple 1/N argument (see 
Fig. 1). We approximated the spatial dependence as a 5- function, though it should be really thought 
of as a function that has support within a distance of order Hubble. For the purposes of proving the 
constancy of at long wavelength, this difference is irrelevant. 

We made the 5-function a scalar function of the local coordinates. This is why there is a factor 
of the determinant of the local metric Ql.iiv^ which is given by (17). The reason why we made the 
5-function a scalar can be understood by noticing that there is a residual gauge-transformation that 
keeps us in the local homogeneous anisotropic coordinates. This is given by a time-independent 
rescaling of the coordinates 

x' M'jx^ . (26) 

This change of coordinates does not affect and A^, but does change and A^, and it can 
be checked that all these variables can be changed independently. In particular, for the choice of 

^This is quite intuitive, but let us give a detailed explanation. The fact that they become uncorrelated in 
time can be quickly understood by realizing that one can rotate the contour of integration in such a way that 
all modes inside the horizon are exponentially damped as e^"^ . Given that the last vertex has to be at about 
horizon crossing, the earlier one cannot be further in the past than about one Hubble time. The fact that 
short modes ks with physical wavenumber greater or equal than Hubble do not contribute to the correlation 
at long physical distances rL{t) ^ can be seen by a simple Fourier space argument: the contribution of 
an e- folding of modes to the correlation at a physical distance rL{t) scales as sin(rL ks/a)/{rL ks/a). Since 
ks /a ^ this contribution vanish as we take rL{t) ^ H~^. 

^An explicit regularization and renormalization using the operator structure of effective field theory of 
single clock inflation was made in [1] and [3] , where the symmetry constraints provided by the effective theory 
turned out to be extremely useful to simplify the calculations. 



= ^^'^ ,7 ""'^ Fa (t, ...)+Fb (t, . . . 
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coordinates in which the gradient of is on the z axis, so that only Al,5 is non- vanishing, we have 
the two independent transformations: 



Cl ^ Cl - « 



(27) 



.^5/3 



V 








2^/3 









-2/3/3 



The right hand side of (25) must therefore be invariant under the residual gauge-transformation. 
This implies that every a::-dependence must be accompanied with a contraction with the local metric. 
The factor of in the local metric must be there because we can choose to re-absorb a isotropic 
rescaling of the coordinates with a rescaling of a. Additionally, this implies that Fa cannot depend 

on and A^, but only on 6Nl^ , all of these vanishing as some power of k\o^/(aH). We 

stress that at diagrammatic level the absence of in the expectation value of the short modes 
appears in general as the result of a cancellation among many different diagrams. The complex 
structure and detailed physical reasons of how this cancellation occurs was studied in [3] . 

The second term represents the contribution from diagrams where xi and X2 are connected by 
low momentum lines (see Fig. 2). We call these Long-Connected (LC) diagrams. At 1-loop in the 2- 
point function calculation, they reduce to the a cut-in-the-side diagrams discussed in [1, 3]. Fb must 
therefore vanish if any of the long wavelength fluctuations vanish. The same gauge transformation 
as in (26) implies that Fb cannot depend on (l and but only on its derivatives, which vanish as 
powers of k\o^/{aH). 

We have taken the dependence of the functions Fa and Fb on (l^ SNlj ... as local in time. 
This is only an approximation. The iterative argument implies that they decay in time with a time 
scale of order Hubble. However, it should be noticed that, because of the ie prescription, the time 
integrals have support of order one Hubble time, and, for the purposes of proving the constancy 
of Cl, we can approximate the dependence as local in time. Fa and Fb can also depend on time, 
but in a way that is proportional to the slow-roll parameters for example through the slow changes 
in H. Also note that these terms include connections through long modes that 'are generated' in 
other vertices that might be outside the volumes around xi and X2 and are connected to both via 
long modes. The appearing inside the Fs account for the entire which is constant over each 
V region (see Fig. 3). 

Let us know take the average over the long modes: 

S^^\xi 



(4^^'^)(xi)4^^'^)(x2))short&long 



X2 



e-3^^F^(t,CL,...)) 



+ 



{FBl{tl)FBl{t2){L{xi)L{x2)) + 



(28) 



where Fbi = OFb/OC^l ^. Let us consider the expectation value in the first line. Since the lower 

order C,l vanish, this is dominated the the expectation value of the determinant of the metric. We 
have 

9, 



(29) 



where we have expanded the argument of the expectation value to the leading order. Since C,l is 
made of lower order modes with wavenumber up to k\Q^^ all of which are constant, this expectation 
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Figure 1: Example of a SC (short-connected) diagram. Red-dashed lines ended by a red-cross 
represent free modes, and two crosses encircled by a blue-dashed ellipse represent a correlation func- 
tion between those two free modes. The dashed-filled ellipse represents possible additional correlation 
functions. This is a SC diagram as, among the various correlation functions present in the diagram, 
there are two involving two short modes at different locations xi and X2 that are contracted among 
each other. If the distance between xi and X2 is much longer than the wavelength and frequency of 
the two modes, the contribution becomes vanishingly small. The points xi and X2 have therefore to 
be very close to each other. This diagram contributes to the first term of (25). 



value is non-perturbatively time- (and space-) independent. The neglected terms decay in time as 
powers of k\o^^/{aH). Analogously the expectation value on the second line decays in time as it is 
made up of lower order C,l that are assumed to be decaying as some power of k\o^/{aH) by the 
iterative argument. The same is true for the higher derivative terms. The details of how the decay 
happens might be dependent on which specific vertex we are considering. In particular, we expect 
generically the suppression to go as {k\ow/ (clH))'^ ^ but note that one factor of {k\ow/{ciH)) suffices 
to make the correlation time independent. The dependence in the first line originates for 
example from short-connected diagrams that are connected additionally by long modes (see Fig. 3). 
Additionally, notice that if there are short vertices around both xi and X2 which are being affected 
by the in the correlation (as opposed of the the line being connected directly to the external (l) 
one should be paying a price in each vertex. 

Let us now compute the two point function of (l- We have 



(a(xi)a(x2)) 



rt rt 

/ dti / dt2 e 



3H{t-ti)^-3H(t-t2) 



(30) 



X2 



a{ti 



^low 

a{ti)H 



FAlitl) + 



+ 



FBl{tl)FBl{t2) 



d 



XI— X2 

aH 



C{xi - X2) -\ 



Mm[ti,t2] 



where C{xi — X2) is the lower-order correlation function of (l, and we have rearranged the terms 
according to their dxi-x2/ i^H) suppression, naming the coefficients of each term as F^..., F^..., .... 
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Figure 2: Example of a LC (long- connected) diagram. Green line represent Green's function. The 
dashed-filled ellipse contains possible additional correlations that connect xi and X2 only though long 
modes. This is a LC diagram as xi and X2 are connected only by low momentum lines. Notice 
that xs must be very close to xi because it is connected to xi by a correlation of short modes. This 
diagram contributes to the second term of (25). 



We have used the slow roll approximation for the Green's function of GL{t, t') ^ Exp(— 3i7(t — t')), 
the subleading corrections being irrelevant. Notice that since C{xi — X2) is approximately scale 
invariant, its spatial derivative is peaked at short distances. 

The time integrals clearly go to zero at late times ^. Notice that the spatial dependence of the 
correlation function consists of the sum of a 5-function and a term that is peaked at short distances 
approximately as 

The same behavior is true when we consider the contribution of n-point correlations involving 
C,L and/or Xl^u- This shows that the correlation functions of C,l and A^,^^ at a given order vanish 
once we take the lower order correlation functions to be time-independent. This proves by iteration 
that C,L and Xl^u have vanishing correlation functions at all orders at late times. They vanish as 
operators. 

Finally, let us consider the long wavelength defined in (20). As given by (12) and (16) in 
Sec. 2, the relationship between C,l and is such that once C,l and A^ vanish so does (^l, implying 
that the correlation functions of C,l vanish to all orders for k\o^^/(aH) <C 1. In particular, since for 
k < fciow, CL{k) = C(^)5 with ({k) being the Fourier transform of the globally defined ( fluctuation, 
the correlation functions of vanish at all orders for k/(aH) <^ 1, as we ultimately wished to show. 

^It might be useful to notice that 

' dtr e-6^(^-^i)e-3^^^ = fl - e'^^'^'-'^-A , 

3 H \ ) 
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Figure 3: Example of a diagram where the dependence of the functions Fa and Fb appearing in 
the average over short modes of {S^^^'^\xi)S^^^'^\x2)) short (top shaded region) is generated by some 
interaction in the past (shaded regions). If the ellipses enclosing xi and X2 contains a correlation 
function made of short modes that connect xi and X2 we obtain a SC diagram, and it contributes 
to the first term of (25). If instead xi and X2 are connected only by long modes, we have a LC 
diagram, and a contribution to the second term of (25). After taking expectation value over the long 
modes, we obtain a contribution respectively to the first and the second line of (28). 



We finally make a comment on the IR divergencies present in the above formula (30). Indeed 
notice that if we try to evaluate (29) perturbatively, we obtain, at leading order, 

d + ^(Cl) + . . . 1 + ^A^ log(fciowi) + . . . , (31) 

where L is the infrared comoving size of the box and where is the tree-level dimensionless power 
spectrum of the C, fluctuations 

(C(fc, t ^ 00)C(fe', t ^ ^))tree-level = {2T^)H^''\k + ^0 ^ , (32) 

where we have neglected the tilt. We see that there is an IR divergence as L ^ oc, hidden in (30) in 
the coeflicients F^..., Fg..., .... Such IR divergencies were found flrst in [14, 15] in the context of the 
long wavelength C, 1-loop power spectrum, but they are present in every ('-correlation function that 
is expressed in comoving coordinates. It was later realized in [2, 16, 17, 18] that they cancel com- 
pletely out once the calculation is performed in local physical coordinates where the long- wavelength 
fluctuations are re-absorbed within the coordinates: for example Vphysicai — Exp(3CL)<^(^)^Komoving- 
Since expressions (28) and (30) are in comoving coordinates, they are also affected by these IR 
divergencies, which are not harmful. More in detail, starting at 2-loops, the IR divergence (31) 
comes from the IR divergence of the short-mode correlation functions at horizon crossing. Notice 
that while for modes longer than the horizon the IR divergence is slow-roll suppressed, this is not 
the case for modes evaluated at horizon crossing. Indeed the three-point function in the squeezed 
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limit is not slow-roll suppressed in this same kinematical regime [3, 11]. The second line of (28) is 
also IR divergent starting at 2-loops as it also involves short mode correlations in the background 
of long wavelength modes. These IR divergencies need to be present in correlation functions in 
comoving coordinates in order to ensure that once the correlation functions are expressed in terms 
of local physical coordinates, every IR dependence on modes much longer than the ones of interest 
is absent [2]. 
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